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MEMORY LSI FAILURE ANALYSIS APPARATUS AND ANALYSIS 

METHOD THEREOF 

BACKGROUND OF THE INVENTION 
The present invention relates to a failure analysis apparatus for 
semiconductor application devices, and in particular relates to a failure 
analysis apparatus for memory LSI. 
5 Description of the Related Art 

Conventional LSI failure analysis devices will now be explained. 
As an example of a memory LSI failure analysis device for improving 
yield and making known the cause of failures, Japanese Patent 
Application Laid-open No. 07-072206 provides an expert system in which 
10 the expertise of a process engineer, a circuit engineer and a layout 
engineer are all incorporated in a personal computer. 

In addition, Japanese Patent Application Laid'Open Nos. 
2000-200814, 2000-321333, 2000 21195, and 2000-187062 provide device 
structures, which distinguish whether or not a defect is a result of the 
15 design by analyzing the types of divisors of the intervals between 
defective bits and their frequencies. 

Moreover, in Japanese Patent Application Laid-open No. 
11-306793, as a defective bits map analysis method, a defect analysis 
method is provided, in which defective bits map is subjected to wavelet 
20 transform, an X-direction/high pass filtered and Y-direction/low pass 
filtered information XH (i, j) is summed for each i in the Y-direction to 
create a histogram, defective addresses X are obtained from those 
addresses i not having their summed value equal to zero, and the 
number of defects is then obtained from the respective absolute values of 
25 those summed values and outputted. 

Since the speed at which the large capacity and the high density 
of memory LSI tend to grow is ever increasing, it has now become 



necessary to respond to the defect analysis of dynamic random access 

memory (DRAM) having 256 megabits, 1 gigabit or more. Added to that 

is the fact that wafer size has also been increasing so that in the future 

they will most certainly have large diameters of 300 mm, and in those 
5 cases, the number of chips and the number of defects to be analyzed 

together with increase drastically. 

In order to solve the above-mentioned problems, high-speed 

techniques have been provided, for example, in Japanese Patent 

Application Laid'Open No. 2000 2 1195, which discloses defect analysis 
io using distributed processing! Japanese Patent Application Laid-Open No. 

2000-187062, which discloses defect analysis using regional divisions! 

and Japanese Patent Application Laid-Open No. 2000-321333, which 

discloses defect analysis algorithms. 

Not only does an increase in the number of defects increase the 
15 load to be analyzed, but it also means that the amount of analyzed 

results output as a result of LSI defect analysis also becomes enormous, 

thus making it difficult to manually judge the cause of the defects. 

For example, considering cases where it is used as a monitoring 

device in order to estimate the reduction in yield from the production line, 
20 since it is impossible to fully confirm the analyzed results manually, the 

defect analysis device itself needs to interpret the analyzed results and 

give off an alarm. 

Accordingly, in future LSI defect analysis devices, not only does 

the analytical ability need to be improved, but also they must have the 
25 ability to automatically interpret the analyzed results. 

There are no conventional LSI defect analysis devices equipped 

with the ability to automatically interpret the analyzed results. 

Furthermore, the device disclosed in Japanese Patent Application 

Laid-Open No. 2000*200814 analyzes the address difference between 
30 defective bits and judges whether or not there is a pattern to the defect 
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distribution; however, the result output by this device is the value of an 
expected value function. In cases where an LSI defect analysis device is 
used as the monitoring device mentioned above, it can be thought of as 
only being useful in setting off an alarm when the rate of patterned 

5 defects exceeds a preset threshold value. Nevertheless, the expected 
value function output from the device disclosed in Japanese Patent 
Application Laid-Open No. 2000-200814 has qualities whereby even if 
the rate of patterned defects remains the same, its value will increase 
according to factor £ therefore, it is not suitable for such usage. 

10 Accordingly, while taking into consideration the above problems 

and limitations, the objective of the present invention is to provide a 
device, method and storage medium, stored with software programs, 
which, when a memory LSI defect analysis device is used as a monitoring 
device to estimate reductions in yield, can shorten the time needed for 

15 full manual interpretation of the obtained results, by automatically 
interpreting the analyzed results obtained, and calculating the period of 
distribution patterns and the mix rate of regular patterned defects. 
Further objectives, characteristics, and benefits of the present invention 
will become apparent to those skilled in the art from the description of 

20 the following preferred embodiments. 

SUMMARY OF THE INVENTION 
The present invention, which achieves above-mentioned objectives, 
comprises: a testing means, which tests LSI; a testing means, which tests 
25 said memory LSI; a data readout means, which reads data that is output 
form said testing means and holds it in memory region of an analytical 
computer; an address difference calculation means, which calculates 
address differences between two defect data; an address difference 
histogram production means, which produces an address difference 
30 histogram based on said address difference! an expected value function 



calculation means, which calculates expected value function T(f) for 
function f of an address difference based on said address difference 
histogram! and a regular patterned defect mix rate calculation 
means, which calculates mix rate of regular patterned defects included in 
a defect distribution from said expected value function. 

The present invention, further comprises a regular patterned 
defect mix rate function calculation means, which finds a regular 
patterned defect mix rate function for every factor f. 

The present invention, further comprises a means of compensating 
a baseline; and a regular patterned defect mix rate function calculation 
means, which finds a regular patterned defect mix rate function based on 
said baseline. 

According to another invention, a memory LSI defect analysis 
method is provided which comprises a first step of testing a memory LSI; 
a second step of reading data associated with defect bits which is 
obtained from results of testing in said first step, and holding the data in 
a memory section of an analytical calculator! a third step of calculating 
address differences between two pieces of defect data; a fourth step of 
producing an address difference histogram based on said address 
differences; a fifth step of calculating expected value function T(f) for 
factor f of an address difference based on said address difference 
histogram; and a sixth step of calculating a regular patterned defect mix 
rate of a defect distribution from said expected value function. 

According to still another invention, a computer program product 
comprising a computer usable medium having computer readable 
program code means embodied in said medium for performing electrical 
testing on memory LSI to be analyzed, by a computer having a memory 
LSI defect analysis apparatus that comprises a memory test system, said 
product comprising! a first computer readable program code means for 
reading data that is obtained as a result of said testing and holding the 



data in memory region of an analytic calculator; a second computer 
readable program code means for calculating address differences between 
two pieces of defect data; a third computer readable program code means 
for producing an address difference histogram based on said address 

5 difference; a fourth computer readable program code means for 
calculating expected value function T(f) for factor f of an address 
difference based on said address difference histogram; and a fifth 
computer readable program code means for calculating mix rate of 
regular patterned defects included in a defect distribution from said 

10 expected value function. 

BRIEF DESCRIPTION OF THE DRAWINGS 
The above-mentioned and other objects, features and advantages 
of this invention will become more apparent by referencing the following 
15 detailed description of the invention taken in conjunction with the 
accompanying drawings, wherein: 

Fig. 1 is a block diagram illustrating the constitution of the first 
embodiment of the present invention; 

Fig. 2 is a block diagram illustrating the constitution of the second 
20 embodiment of the present invention; 

Fig. 3 is a block diagram illustrating the constitution of the third 
embodiment of the present invention; 

Fig. 4 is a flowchart illustrating the constitution of the first 
embodiment of the present invention; 
25 Fig. 5 is a flowchart showing the process flow of regular patterned 

defect mix rate calculation means 16 according to the first embodiment of 
the present invention; 

Fig. 6 is a flowchart showing the process flow of regular patterned 
defect mix rate function calculation means 21 according to the second 
30 embodiment of the present invention; 




Fig. 7 is a flowchart showing the process flow of regular patterned 
defect mix rate calculation means 32 according to the third embodiment 
of the present invention! 

Fig. 8 is a diagram for describing the first embodiment of the 
5 present invention, and is a diagram showing an example of defect data 
including regular patterned defects at period 10; 

Fig. 9 shows a graph of expected value function T(£) for defect data 
A in Fig. 8; 

Fig. 10 illustrates a graph of expected value function T(f) for 
10 defect data B in Fig. 8; 

Fig. 11 is a graph of expected value function T(f) for defect data C 
in Fig. 8; 

Fig. 12 shows a graph of expected value function T(£) for defect 
data D in Fig. 8; 

15 Fig. 13 is a graph of expected value function T(f) for defect data E 

in Fig. 8; 

Fig. 14 is a diagram showing a graph in which the relationship 
between the ratio of (T ma x - l)/(fmax - l) and the ratio of the number of 
regular patterned defects over the total number of defects is shown! 
20 Fig. 15 is a diagram showing a graph in which the relationship 

between the ratio of (Tmax - l)/(fmax - l) and the square of the ratio of the 
number of regular patterned defects over the total number of defects is 
shown! 

Figs. 16A and 16B are a diagram describing baseline 
25 compensation! 

Fig. 17 shows an example of a graph of expected value function 
T(£); and 

Fig. 18 is a diagram showing a graph of regular patterned defect 
mix rate function MR(f) corresponding to the defect data used in Fig. 17. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
The preferred embodiments of the present invention will now be 
described. In a first embodiment of the present invention, which will be 
described in detail later, a memory LSI defect analysis device, which 

5 comprises a memory test system that performs electrical testing of 
memory LSI, is comprised of testing means 11 (see Fig. l), which tests 
memory LSI; data readout means 12, which reads out the data of defect 
bits that have been read out from the testing means and records it; 
address difference calculation means 13, which calculates the iaddress 

10 difference between two defect data; address difference histogram 
preparation means 14, which creates a histogram of address differences 
based on the address differences; expected value function calculation 
means 15, which calculates expected value function T(£) for factor f based 
on the address difference histogram; and, regular patterned defect mix 

15 rate calculation means 16, which calculates the mix rate of regular 
patterned defects that are included in the defect distribution by the 
expected value function. 

In more detail, regular patterned defect mix rate calculation 
means 16, as shown in Fig. 5, 

20 (a) finds maximum value T ma x of expected value function T(£) for 



factor f, 

(b) finds value fmax of fat the time of maximum value T ma x, 

(c) decides if maximum value Tmax is larger than 1, 

(d) decides if larger than 1 then "regular patterned defect", or if 1 
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or smaller then "irregular patterned defect", and 
(e) when regular patterned defect, calculates mix rate with the 



following formula- 



regular patterned defect mix rate = 




In another embodiment of the present invention, instead of 



regular patterned defect mix rate calculation means 16, regular 
patterned defect mix rate function calculation means 21 may be provided, 
which calculates regular patterned defect mix rate function MR(f) from 
expected value function Tit). 

As shown in Fig. 6, regular patterned defect mix rate function 
calculation means 21 

(a) selects factor f and finds the value of expected value function 
T(£) for the factor f, 

(b) decides if T(£) is larger than 1 or not, 

(c) if T(£) is larger than 1, then decides the defect distribution has 
a regular pattern at period f and finds the value of regular 
patterned defect mix rate function MR(£> using the following 
formula : 



(d) if T(f) is 1 or less, then decides the defect distribution does not 
have a regular pattern at period f and makes 

regular patterned defect mix rate function MR(f) = 0 , and 

(£) confirms whether or not regular patterned defect mix rate 
function MR(f) has been found for all factors f, and if finished 
then ends one string of steps, if not returns to step (a). 

In another embodiment of the present invention, baseline 
compensation means 31, which allows calculation of a more accurate mix 
rate, and regular patterned defect mix rate function calculation means 32, 
which calculates a regular patterned defect mix rate function 
corresponding to a baseline compensated by the baseline compensation 
means 31 are provided. 

As shown in Fig. 7, regular patterned defect mix rate function 
calculation means 32 

(a) finds total number of defects n using defect calculation 



MR(£> 
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processing, 

(b) selects factor f, 

(c) finds the value of expected value function T(f) for that f, and 
decides whether or not the value is (n-f)/(n-l) or greater and 

5 that the number of defects n is equal to or greater than factor f, 

(d) when the conditions of T(f) > (n - f)/(n - l) and n > f have been 
met, decides defect distribution has a regular pattern at period 
f, and finds regular patterned defect mix rate function MR(f) 
using the following formula 

MR(fi = / <"-'PX/)- + / , 

(e) when the above-mentioned conditions are not met, decides 
defect distribution does not have a regular pattern at period f, 
and makes 

regular patterned defect mix rate function MR(£) = 0, 
15 (£) confirms whether or not regular patterned defect mix rate 

function MR(£) has been found for all of the factors f, and if 
finished then ends one string of steps, if not returns to step (b). 
In the present invention, data readout means 12, address 
difference calculation means 13, histogram preparation means 14, 
20 expected value function calculation means 15, regular patterned defect 
mix rate function calculation means 16, regular patterned defect mix rate 
function calculation means 21, baseline compensation means 31, and 
regular patterned defect mix rate function calculation means 32 are able 
to execute their processes and functions using the program run in the 
25 computer which comprises the memory defect analysis device. In this 
case, in the computer that contains a storage medium readout device, the 
storage medium (e.g. floppy disk (FD), hard disk drive (HDD), magnetic 
tape (MT), compact disk ROM (CD-ROM), digital versatile disk (DVD), 
semiconductor memory, etc.), which records the program, reads it out to 




the computer via the storage medium readout device. Alternatively, the 
program can be downloaded to the computer via a communication means, 
then be loaded into the main memory and be executed in order to put 
into operation the memory defect analysis device of the present 
5 invention. 

The embodiments of the present invention will now be described 
below while referencing the Figures. 
First embodiment: 

Fig. 1 is a diagram showing the structure of a first embodiment of 
10 the memory LSI defect analysis device of the present invention. Fig. 4 
is a flowchart showing the process flow according to the first embodiment 
of the present invention. 

As shown in Figs. 1 and 4, in the first embodiment of the present 
invention, testing means 11 performs electrical testing on the memory 
15 LSI to be tested, the results are then output as bitmap data to a memory 
device. The bitmap data is read out by the data readout means 12, and 
the coordinate data of each defect bit is held in the memory of analytical 
computer (Step 41). 

In the address difference calculation means 13, two defect bits, a 
20 and b, of the defect bits that have been read are selected, and the 
difference of their addresses d(a, b) is calculated (Step 42). The method 
of calculating this address difference varies depending on the type of 
analysis; however, when X address analysis is performed, the address 
difference can be shown as the absolute value of the difference of the X 
25 coordinates of defect bits a and b. In other words, if the X coordinates of 
a and b are given as a(x) and b(x), respectively, then 

d(a,b)= | a(x)-b(x) | ...(1). 
On the other hand, when Y address analysis is performed, the 
address difference is calculated to be the absolute value of the difference 
30 of the Y coordinates of defect bits a and b. That is, if the Y coordinates 



of a and b are given as a(y) and b(y), respectively, then 

d(a,b)= I a(y)-b(y) I ...(2). 

In the histogram preparation means 14, summation is performed 

in order to prepare histogram H(d) for address difference d. Specifically, 
5 for defect bits a and b, H(d(a, b)) is incremented by 1 (Step 43). 

A completed address difference histogram H(d) can be obtained by 

performing the operations of address difference calculation means 13 and 

histogram preparation means 14 for each respective pair of defect bits 

(Steps 42 to 44). 

10 In the expected value function calculation means 15, the 

calculation of expected value function T(f) is performed based on the 
address difference histogram H(d) (Step 45). 

It should be noted that the expected value function T(i) has been 
defined as equation (3) shown below (reference: Japanese Patent 

15 Application Laid-Open No. 2000 200814). 

TOD = y£w(/) ...(3) 
(N-ux) 

where, 

Zm(/) : the number of combinations of defect bits where address 
difference has / as a factor, 
20 N ' total number of combinations of defect bits, and 

ux- the number of combinations of defect bits where address difference 
becomes 0. 

When max(d) denotes the maximum address difference of defect 
bits, n denotes the number of defect bits, and j denotes the counter, the 
25 following equations are given: 

N = ^2),and 
2 

ux = H(0) 

Therefore, expected value function T(f) can be found from address 
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difference histogram H(d) using equation (4) shown below. 

T(£> = 2tf{n(n-l)-2H(0)}x2(fj) ...(4) 

Here, the summing of 2 H(fj) is performed for j where j = 1 to fj ^ 

max(d). 

5 In the regular patterned defect mix rate function calculation means 

16, the mix rate of regular patterned defects is calculated based on 
expected value function T(£> (Step 46). It is possible to find the 
relationship between the expected value function T(f) needed for this, 
and the mix ratio of regular patterned defects in the following manner' 

10 The mix ratio of regular patterned defects is defined as the 'ratio of 

the regular patterned defects to the total number of defects'; however, 
there exists the following relationship between the expected value 
function T(f) for defective data, which includes regular patterned defects 
with period X , and the mix rate of regular patterned defects. 

15 Figs. 9 through 13 show the results from calculating the expected 

value function T(f) for defect data A through E, respectively, which are 
shown in Fig. 8 and include regular patterned defects with period 10. 

As shown in Fig. 8, data A has a defect distribution where each 
coordinate of address coordinates 0 to 1000 is a defect bit, and from 1000 

20 to 50000, there is a defect bit every tenth coordinate (e.g., 1010, 1020, ... , 
49990, 50000). The total number of defect bits is 5900, and out of that, 
4900 are regular patterned defects with period 10, therefore, the regular 
patterned defect mix rate is 
4900/5900 = 83.05%. 

25 Data B has a defect distribution where each coordinate of address 

coordinates 0 to 10000 is a defect bit, and from 10000 to 50000, there is a 
defect bit every tenth coordinate (e.g. 10010, 10020, ... , 49990, 50000). 
The total number of defect bits is 14000, and out of that, 4000 are regular 
patterned defects with period 10, therefore, the regular patterned defect 

30 mix rate is 




4000/14000 = 28.57%. 

As for Data C, it has a defect distribution where each coordinate of 
address coordinates 0 to 20000 is a defect bit, and from 20000 to 50000, 
there is a defect bit every tenth coordinate (e.g., 20010, 20020, ... , 49990, 
5 50000). The total number of defect bits is 23000, and 3000 out of 23000 
are regular patterned defects with period 10. Therefore, the regular 
patterned defect mix rate is 
3000/23000 = 13.04%. 

Data D has a defect distribution where each coordinate of address 
10 coordinates 0 to 30000 is a defect bit, and from 30000 to 50000, there is a 
defect bit every tenth coordinate (e.g., 30010, 30020, ... , 49990, 50000). 
The total number of defect bits is 32000, and out of that, 2000 are regular 
patterned defects with period 10, therefore, the regular patterned defect 
mix rate is 
15 2000/32000 = 6.25%. 

Data E has a defect distribution where each coordinate of address 
coordinates 0 to 40000 is a defect bit, and there is a defect bit every tenth 
coordinate (e.g., 40010, 40020, ... , 49990, 50000) from 40000 to 50000. 
The total number of defect bits is 41000, and out of that, 1000 are regular 
20 patterned defects with period 10. The regular patterned defect mix rate 
is therefore 

1000/41000 = 2.44%. 

At the same time, as shown in the graphs of the respective 
expected value function T(£) in Figs. 9 through 13, the height of the 
25 respective peaks gradually decreases as follows : 

in Fig. 9, which is the analyzed result of data A, T(10) = 7.207, 
in Fig. 10, which is the analyzed result of data B, T(10) = 1.734, 
in Fig. 11, which is the analyzed result of data C, T(10) = 1=153, 
in Fig. 12, which is the analyzed result of data D, T(l0) = 1.035, 

30 and 




in Fig. 13, which is the analyzed result of data E, T(10) = 1.005. 
Now, the expected value function T(f) has the following properties. 

(1) When the defect distribution has only regular patterned defects 
with period X , T(f) has its peak at f = (k = 1, 2, 3, ...), and its 
maximum value is T( X ) = X , when f = X . 

(2) When the defect distribution is completely irregular patterned 
defects, T(f) = 1 for all f. 

It can be regarded that in (l) above, the regular patterned defect 
mix rate is 100%, and in (2) above, the regular patterned defect mix rate 
is 0%. Accordingly, the mix rate of regular patterned defects with period 
fmax can be expressed as the value of 

(T max 

-l)/(£nax-D, 

where Tn, ax is the maximum value of expected value function T(f), and 
fmax is the value of factor f when expected value function T(f) has the 
maximum value T max* 

In the examples of data A through data E, since fmax = 10, the 
relationship of the value of (T(l0)-l)/(l0-l) of each data and the value of 
the ratio of the number of regular patterned defects over the total 
number of defects is as shown in the graph in Fig. 14. The graph in Fig. 
14 shows the relationship between (T ma x-l)/(f«nax-l) and the ratio of the 
number of regular patterned defects over the total number of defects. 

Note that the data that has been plotted also utilizes 35 pieces of 
data differing from the respective mix rates in addition to data A through 
E. 

As it is clear from the graph in Fig. 14, it can be understood that 
there exists the relationship of a second order equation between the 
value of (T(10)-1)/(10-1) in each data and the value of the number of 
regular patterned defects over the total number of defects. 

Fig. 15 is the graph where the relationship between the value of 
(T(10)-1)/(10-1) for each data and the squared value of the ratio of the 
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number of regular patterned defects over the total number of defects is 
plotted. According to the graph in Fig. 15, it can be understood that 
(T(10)-1)/(10-1) 

= {(the number of regular patterned defects)/(the total number of 
defects)) 2 -..(5) 

Accordingly, using maximum value T max of expected value function 
T(£>, and value fmax of fat that time, the following formula is obtained: 



regular patterned defect mix ra te = JLS^zl ...( 6 ) 

y / max- 1 

From the above, the process of regular patterned defect mix rate 
10 calculation means 16 can be explained by the flowchart in Fig. 5. 

Specifically, the maximum value T max of expected value function 
T(f) is first found in Step 51, then fmax, or the value of factor f when the 
expected value function T(f) takes its maximum value Tmax is obtained in 
Step 52. 

15 Next, whether or not T max is greater than one is decided in Step 

53; if it is greater than 1, then 'regular pattern distribution (Step 54); if 
it is not, then 'irregular pattern distribution' (Step 55). 

When it is regular pattern distribution, the mix rate is found in 
Step 56 using the following formula: 



20 regular patterned defect mix rate = J ^ max 1 ...(7). 

V / max - 1 

Second embodiment: 

The second embodiment of the memory LSI defect analysis device 
according to the present invention will be described. Fig. 2 is a diagram 
showing the structure of the second embodiment of the memory LSI 
25 defect analysis device according to the present invention. Fig. 7 is a 
flowchart showing the process flow of the second embodiment of the 
present invention. 

As shown in Fig. 2, the memory LSI defect analysis device of the 




second embodiment comprises storage medium 11, data readout means 
12, address difference calculation means 13, histogram production means 
14, and expected value function calculation means 15 as with the first 
embodiment described above as shown in Fig. 1, and in addition to these, 
5 regular patterned defect mix rate function calculation means 21 is 
provided. 

The second embodiment of the present invention calculates 
regular patterned defect mix rate function MR(f) from expected value 
function T(f) in the regular patterned defect mix rate function calculation 
10 means 21. 

In the first embodiment described above, the formula that 
calculates the regular patterned defect mix rate can be applied not only 
for f when T(f) is maximum, but for all values of f. Therefore, regular 
patterned defect mix rate function MR(£)can be thought of as a function 
15 of f as with T(f), and it can be defined by the following equation: 

M R(f) - (8) 

As a result, the process of regular patterned defect mix rate 
function calculation means 21 can be explained using the flowchart 
shown in Fig. 6. 

20 In the regular patterned defect mix rate function calculation 

means 21, factor f is first selected in Step 61, the value of expected value 
function T(D for that f is then found, and it is decided in Step 62 whether 
or not that value is greater than 1. If it is greater than 1, then the 
defect distribution is determined to be regularly patterned at period f 

25 (Step 63), and regular patterned defect mix rate function MR(f) is found 
in Step 64 using the following formula: 

On the other hand, if the value of T(£) is equal to or less than 1, 




then the defect distribution is determined, in Step 65, not to be regularly 
patterned at period f, and regular patterned defect mix rate function 
MR(£) becomes equal to 0 (Step 66). 

Finally, it is confirmed in Step 67 that regular patterned defect 
5 mix rate MR(f) has been found for every f, and if finished then one string 
of steps ends, if not then it returns to Step 61. 
Third embodiment: 

The third embodiment of the memory LSI defect analysis device 
according to the present invention will be described. Fig. 3 is a diagram 
10 showing the structure of the third embodiment of the memory LSI defect 
analysis device according to the present invention. Fig. 7 is a flowchart 
showing the process flow of the third embodiment of the present 
invention. 

As shown in Fig. 3, the memory LSI defect analysis device of the 
15 third embodiment comprises storage medium 11, data readout means 12, 
address difference calculation means 13, histogram production means 14, 
and expected value function calculation means 15 as with the first 
embodiment described above as shown in Fig. 1, and in addition to these, 
baseline compensation means 31 and regular patterned defect mix rate 
20 function calculation means 32 are provided. 

In the third embodiment of the present invention, the baseline 
compensation means 31 performs baseline compensation during 
calculation of regular patterned defect mix rate function MR(f). The 
regular patterned defect mix rate function calculation means 32 then 
25 calculates, based on that compensation, regular patterned defect mix 
rate function MR(£). 

First, the baseline compensation means 31 will be described. 
With the equation for the regular patterned defect mix rate in the 
first and second embodiments described above, the expected value 
30 function T(£> of irregularly distributed defects becomes equal to 1 when 




there is an unlimited number of defects; or in other words, it is assumed 
that T(f) becomes equal to 1 for all factors f. 

However, since there is inevitably a limit to the number of defect 
data in the actual object that is analyzed, in order to calculate a more 
5 accurate mix rate, baseline compensation in the form of T(f) = 1 must be 
performed. 

In order to perform this, the value of the expected value function 
for irregularly distributed defects in the case where there was a limit to 
the number of defects, is first found, and then the mix rate formula may 
10 be defined in accordance with it. 

In the following, the value of the expected value function for 
irregular patterned defects having a limited number of defects is found. 

To begin with, it is first assumed to have an irregular distribution 
of n-number of defects where a defect bit exists at each coordinate within 
15 coordinates 1 to n. 

Accordingly, 

there are n-2 defect pairs having a distance of 2- (l, 3), (2, 4), ... , 
(n-2, n), and 

there are n-3 defect pairs having a distance of 3: (l, 4), (2, 5), ... , 
20 (n-3, n), 

and similarly, 

there are n~d defect pairs having a distance of d: (l, d+l), (2, d+2), 
... , (n-d, n), and 

there are n-fi defect pairs having a distance of fi:(l, fi+l), (2, fi+2), 
25 ... , (n-fi, n). 

At this point, if the number of defect pairs having multiples of 
distance f are counted, the counted value becomes the sum of the number 
of pairs in terms of distances: f, 2f, 3f, .... If it is defined as k = [n/f] ( the 
largest integer not exceeding n/f), then the following equation can be 
30 given: 
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the number of pairs each having a distance that is a multiple of f 

_ k(2l1-fk-f) (g.^ 

2 

On the other hand, the total number of defect pairs is equal to the 
5 number of combinations n C2 = n(n-l)/2. Accordingly, regular patterned 
defect ratio P(f) is given by, 

«(« - 1) 2 
_ k(2n-flc-f) (9 . 2) 

Therefore, expected value function T(f) becomes 
10 T(f) = fP(f) = ^( 2// --^-/) ...(10). 

- 1) 

At this point, since, if n is sufficiently large, it is regarded as 
n = fk 
obtaining 

T(£> = = ...(11). 

n(p - 1) w - 1 

15 As a result, it is possible to compensate the baseline from what is 

shown in Fig. 16A (before compensation, T(f) = l) into what is shown in 
Fig. 16B. Fig. 16B shows expected value function T(f) once it has been 
compensated using equation (ll) above, based on the relationship to the 
number of defects n for factor f. 

20 In the regular patterned defect mix rate function calculation 

means 32, regular patterned defect mix rate function MR(f) is calculated 
in the manner described below from expected value function T(i) based 
on the above baseline. 



10 



20 



MR(£>2= = (n-l)n/)-n + f (u) 
y n-f »(/-!) 
n- 1 

Therefore, 



However, when T(f) < " ^ and n < f , then 

MR(£>=0 ...(14). 
Accordingly, the processes of regular patterned defect mix rate 
function calculation means 32 can be described using the flowchart 
shown in Fig. 7. 

Namely, the total number of defects is found in Step 71, the factor 
f is then selected in Step 61. The value of expected value function T(£) 
for that f is found, and it is decided whether or not that value is equal to 
or greater than (n-£)/(n-l) as well as whether or not the number of defects 
n is equal to or greater than factor fin Step 72. 

If the conditions of Step 72 have been met, then it is decided that 
15 the defect distribution is regularly patterned at period f in Step 63, and 
then the value of regular patterned defect mix rate function MR(£) is 
found in Step 73 using the following equation: 



M R(0= + / 



On the other hand, if the conditions of Step 72 have not been met, 
20 then it is decided in Step 64 that the defect distribution is not regularly 
patterned at period f, and thereby making regular patterned defect mix 
rate function MR(£) = 0 in Step 66. 

Next, and finally, it is confirmed in Step 67 whether or not regular 
patterned defect mix rate function MR(f) has been found for every f, and 
25 if finished then one string of steps ends. If not, then it returns to Step 




61. 

Expected value function T(f) has a tendency to increase with 
factor f, even if the mix rate of regular patterned defects does not change. 
However, with regular patterned defect mix rate function MR(f), since it 
5 does not have such properties, it can also be put into use as a new index 
for estimating the cause of defects. 

For example, in Fig. 17, expected value factor T(£) is found up to f 
= 128 for certain defect data, and that data is shown as a graph. 
Looking at Fig. 17, it is noticeable that, since expected value function T(£) 
10 is at its maximum value 3.7 at the point where f = 114, at first glance it 
seems as if there are the greatest number of regular patterned defects at 
period 114. 

In contrast, in Fig. 18, regular patterned defect mix rate function 
MR(f) is found for the same data as one used above, and shown as a 
15 graph. Looking at Fig. 18, it is clear that there are regular patterned 
defects at periods (f = 14, 19, 31, 38) having higher mix rates than f = 
114. 

As described above, in the conventional LSI defect analysis device, 
assuming it is used during the LSI development stage or trial 

20 manufacturing stage, it only outputs expected value function T(£> of 
regular patterned defects. The memory LSI defect analysis device 
according to the present invention is able to calculate regular patterned 
defect mix rate function MR(£> from the expected value function T(f) of 
regular patterned defects. Therefore, it is possible to operate the LSI 

25 defect analysis device in cases, for example where it is used in a 
semiconductor production line as a monitoring device which gives off an 
alarm when the ratio of regular patterned defects exceeds a certain 
preset threshold value, making it advantageous. 

Furthermore, the value of expected value function T(f) based on 

30 factor f tends toincrease as factor f becomes larger even if the mix rate of 




regular patterned defects does not change. However, since expected 
value defect mix rate function MR(£) does not have such property, with 
the present invention, it is easier to spot periods having a high mix rate 
of regular patterned defects, it can also be put into use as a new index for 
estimating the cause of defects. 



